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Singlet fission can generate an exchange-coupled quintet triplet pair state "TT, which could lead to the realization 
of quantum computing and quantum sensing using entangled multiple qubits even at room temperature. However, 
the observation of the quantum coherence of "TT has been limited to cryogenic temperatures, and the fundamen- 


Copyright © 2024 tie 
Authors, some rights 
reserved; exclusive 
licensee American 
Association for the 
Advancement of 
Fujiwara’, Nobuo Kimizuka’”, Science: No claim to 
original U.S. 
Government Works. 
Distributed under a 
Creative Commons 
Attribution 
NonCommercial 
License 4.0 (CC BY-NC). 


tal question is what kind of material design will enable its room-temperature quantum coherence. Here, we show 
that the quantum coherence of singlet fission-derived *TT in a chromophore-integrated metal-organic framework 
can be over hundred nanoseconds at room temperature. The suppressed motion of the chromophores in ordered 
domains within the metal-organic framework leads to the enough fluctuation of the exchange interaction neces- 
sary for "TT generation but, at the same time, does not cause severe °TT decoherence. Furthermore, the phase and 
amplitude of quantum beating depend on the molecular motion, opening the way to room-temperature molecu- 


lar quantum computing based on multiple quantum gate control. 


INTRODUCTION 

We are currently experiencing a second quantum revolution. Quan- 
tum information science (QIS) uses the principles of quantum 
mechanics for various applications such as quantum computing (1, 2), 
quantum communication (3, 4), and quantum sensing (5-7). The 
most basic information unit in QIS is a quantum bit (qubit). Among 
various types of qubit, those using molecular materials generally 
have the following advantages: the ability to precisely create a spe- 
cific qubit structure, the ability to control qubit properties by chang- 
ing the chemical structure, and the scalability to integrate many 
qubits (8-10). While most qubits require cryogenic temperatures, 
some molecular qubits can be initialized by photoexcitation (11, 12), 
making them particularly promising for quantum computing at 
room temperature and quantum sensing of living systems. 

Among the photogenerated molecular qubits, singlet fission (SF) 
has the unique feature of being able to generate high-spin quintet state, 
which is useful for advanced logic gate operations and entanglement- 
enhanced sensing (13-26). SF is a process in which two triplet exci- 
tons (Tı) are generated from one photoexcited singlet (S1) and has 
been studied mainly in the energy field to improve the efficiency of 
solar cells (27-29). On the other hand, many studies on the mecha- 
nism of SF have revealed that the triplet pair with singlet multi- 
plicity TT, in which two T; are strongly coupled by an exchange 
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interaction (J), arises from Sı, and the following spin evolution of 
the triplet pair leads to the formation of a four-spin entangled quin- 
tet multiexciton "TT (13-26). SF-derived °TT can play an important 
role in the future development of molecular quantum computing 
using multiple qubits such as CNOT gate operation. “TT can not 
only be generated by light but also be detected by light through the 
inverse process of SF, triplet-triplet annihilation, and thus has the 
potential to be addressed in a single spin level (15). 

To maximize the potential of these SF-derived °TT as qubits, it is 
crucial to generate quantum spin coherence in the quintet sublevels 
by microwave manipulation at room temperature. However, only a 
few cases of “TT quantum coherence have been observed so far, and 
these cases are limited to cryogenic temperatures of 75 K or lower 
(13-19). This poses a major challenge for the development of SF- 
derived "TT qubits. Most of the reported examples are based on 
molecular crystals of tetracene derivatives, where the two triplet 
excitons move apart, weakening the exchange interaction J and 
causing mixing between 'TT and °TTT, and then the triplets come 
close to each other again, producing a strongly exchange-coupled 
°TT. In this mechanism, delicate control of molecular packing is 
required to make the interaction between chromophores strong 
enough to cause SF but weak enough to suppress the triplet diffu- 
sion (15). In particular, at room temperature, the triplet diffusion is so 
pronounced that the °IT pairs quickly dissociate to form noncor- 
related T + T states, and no coherence can be observed (13). 

Another mechanism for SF-based "TT generation is the modula- 
tion of exchange interactions between chromophores by molecular 
motion, causing the conversion of ITT to °TT (19-22). Such a con- 
formational change has been observed in intramolecular SF systems 
consisting of discrete molecules of covalently linked acene units 
dispersed in a low-temperature glass matrix and may also occur in 
acene dimers doped in molecular crystals (18-20). However, while 
molecular motion is necessary for the generation of "TT, decoherence 
may be induced by fluctuations in the zero-field splitting interaction 
(19). That is, it remains unclear what kind of molecular motion 
can simultaneously achieve efficient “TT generation and noise- 
suppressed qubits. 
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Here, we report the first room-temperature observation of the 
quantum coherence of SF-derived °TT by the suppressed molecular 
motion in a metal-organic framework (MOF). MOFs are nanoporous 
crystalline materials composed of metal ions and organic ligands 
(30-37). By integrating chromophores into the ligand, the distance 
and angle between chromophores can be precisely regulated in 
MOFs (31-33). Unlike dense molecular crystals, MOFs have nanoscale 
voids in the crystals, which allow the chromophores to move, and 
this movement can be controlled by the topology of the network and 
the local molecular density around the chromophores. As a proof of 
concept, we prepare a UiO-type (37) pentacene-based MOF (Pn-MOF), 
which is synthesized by combining diamagnetic Zr ions with a 
dicarboxylate ligand containing pentacene (38-40), the most repre- 
sentative chromophore exhibiting exothermic SE Note that this 
Pn-MOF is the first example of an MOF exhibiting SF to the best of 
our knowledge. The noninterpenetrated UiO-type structure can 
prevent m-stacking between pentacene planes, which gives the pen- 
tacene units enough motion to make the conversion from 'TT to 
STT (Fig. 1). The dense integration in the Pn-MOF makes the pen- 
tacene motion sufficiently suppressive that the quantum coherence 
of °TT generated by microwave irradiation was observed for over 
hundred nanoseconds even at room temperature. Furthermore, the 
phase and amplitude of quantum beating can be dependent on the 
molecular motion. The decay of the echo signal from the quintet 
observed by pulsed electron paramagnetic resonance (EPR) was 
considerably longer than the decay of the quintet-derived signal in 
continuous-wave time-resolved EPR (CW-TREPR) measurements, 
suggesting that two components with different degree of order exist 
and the highly crystalline domains present as minor components in 
the disordered MOF structure exhibit the long quantum coherence. 
These results reveal how molecular motions should be designed to 
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achieve quantum coherence of multilevel distinct quantum systems, 
referred to as qudits, even at room temperature. This provides a 
fundamental guideline for the future development of molecular- 
based QIS. 


RESULTS 

As a dicarboxylate ligand containing a pentacene moiety, we used 
4,4’-(pentacene-6,13-diyl)dibenzoic acid (PDBA; Fig. 2A), which 
we had previously developed for other purposes (36, 41). The energy 
levels of excited states of PDBA were estimated from time-dependent 
density functional theory calculations. B3LYP/6-31G (d,p) and 
@B97X/6-31G (d,p) were used from the ground-state optimization 
and excited-state calculations, respectively. It was confirmed that 
PDBA has the similar Sı (2.34 eV) and T; (0.93 eV) energy levels 
compared with those of pentacene (Sj, 2.45 eV; T1, 0.89 eV) (42), 
and exothermic SF is expected to occur for PDBA as well. We have 
previously confirmed that aggregates of PDBA in water exhibit SF 
(41). Pn-MOF with pentacene integrated in a UiO-type structure 
was synthesized by a hydrothermal method using ZrCl, and PDBA 
in the presence of acetic acid as the modulator in deaerated N,N’- 
dimethylformamide (DMF) at 403 K. The coordination of the 
carboxylic acid of PDBA to the Zr-based cluster was confirmed by 
the shift to lower wave number of the peak originating from car- 
bonyl stretching from 1685 to 1595 cm™' (asymmetric) and 1540 cm~! 
(symmetric) (fig. S1) (43). The structure of Pn-MOF was confirmed 
by powder x-ray diffraction (PXRD) measurements. The peak posi- 
tions of the experimental PXRD pattern of Pn-MOF matches well with 
those of the simulated pattern (Fig. 2B). The simulated structure with 
UiO-68-type topology using PDBA asa ligand is shown in Fig. 2D. The 
center-to-center distance between the nearest pentacene moieties is 
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Fig. 1. Quintet multiexciton generation and its quantum coherence in an MOF. (A 


TT, conformation 
(larger exchange coupling) 


) SF followed by spin evolution via modulation of exchange interactions results in 


the formation of quintet multiexciton °TT in an MOF. (B) Suppressed chromophore motion in an MOF enables not only the spin conversion from 'TT to °TT but also the 


quantum coherence of "TT over hundred nanoseconds even at room temperature. 
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Fig. 2. Characterizations of Pn-MOF. (A) Chemical structure of PDBA. (B) Simulated PXRD pattern of Pn-MOF (black line) and experimental PXRD patterns of Pn-MOF (red 
line), Pn-MOF in paraffin (green line), and paraffin (blue line). a.u., arbitrary units. (C) UV-vis absorption spectra of 1 mM PDBA in DMF (black line), Pn-MOF (red line), neat 
PDBA solid (blue line), and Pn-MOF in paraffin (green line). (D) Simulated structure of Pn-MOF with UiO-68-type topology. (E) The nearest and second nearest pentacene 


moieties in Pn-MOF. 


1 nm, and the dihedral angle is 120°. The distance between the nearest 
carbons of the nearest pentacene pair is 0.35 nm, which is close 
enough for SF to occur. The distance between the second nearest 
pentacene sites is 1.5 nm, and the dihedral angle is 90°. The rela- 
tively broad PXRD peaks of Pn-MOF suggest that the structure is 
partially disordered (44). Scanning electron microscopy images of 
Pn-MOF showed a rounded morphology with a grain size of ~50 to 
100 nm (fig. S2), suggesting not high crystallinity of Pn- MOF, con- 
sistent with the PXRD results. In Pn-MOE, pentacene moieties do 
not form a-stacking with each other, which is expected to suppress 
immediate deactivation of triplet dimer by triplet-triplet annihilation 
and to fluctuate the exchange interaction by changing the dihedral 
angle between pentacene units. These moderate interpentacene 
interactions were also confirmed by ultraviolet (UV)-visible (vis) 
absorption measurements (Fig. 2C). The 0-0 absorption peak of 
the m-x* transition of the pentacene moiety was located at 
602 nm in DMF solution of PDBA, while the absorption peak of 
Pn-MOF showed a slight broadening and red-shift to 612 nm, 
indicating the presence of excitonic interactions between the penta- 
cene moieties. On the other hand, in the neat PDBA solid, the ab- 
sorption peak was located at a longer wavelength of 622 nm with 
more broadening. This indicates that the interaction is properly 
modulated in Pn-MOF by preventing a-stacking between penta- 
cene moieties. 

From thermogravimetric analysis, it was found that Pn-MOF has 
high thermal stability up to about 573 K and also adsorbs water 
from the air into its nanopores (fig. $3). The existence of nano- 
pores in Pn-MOF was also confirmed by nitrogen gas adsorption 
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measurements at 77 K after activation of Pn-MOF at 373 K under 
vacuum (fig. S4). The adsorption isotherm of Pn-MOF showed a steep 
rise in the low-pressure region, which is typical for microporous 
materials and classified as International Union of Pure and Applied 
Chemistry (IUPAC) type I. Gradual rise was also observed in the 
high-pressure region, implying the broad distribution of the pore 
sizes and the existence of defects. Once the nanoporous structure of 
Pn-MOF was confirmed, Pn-MOF was soaked in paraffin. Paraffin 
is nonvolatile, useful for dispersing MOF solids, and is sometimes 
used for time-resolved transient EPR measurements as an inert 
matrix. The PXRD pattern of activated Pn-MOF soaked in paraffin 
retained the original peaks, indicating that the UiO-68-type topology 
was preserved (Fig. 2B). The UV-vis absorption peak shift of 
Pn-MOF was very small (612 to 615 nm) before and after the 
paraffin soaking (Fig. 2C). Therefore, it is likely that paraffin in the 
MOF nanopores does not have any strong interactions with chro- 
mophores and interchromophore interactions were not affected by 
the introduction of guest paraffin molecules. 


SF in Pn-MOF 

To investigate the SF dynamics in Pn-MOF in paraffin, we carried 
out pump-probe-type transient absorption spectroscopy (TAS) in 
the 460- to 550-nm range (Fig. 3), commonly used for tracking SF 
dynamics in pentacene-based systems. Upon excitation by an ultra- 
short pump pulse at 600 nm, we observed a broad transient absorp- 
tion spectrum around 450 to 500 nm, typically attributed to the 
feature of S-S, absorption of pentacene chromophore. The spectrum 
subsequently transformed to a different shape with a peak at 530 nm 
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Fig. 3. TAS measurements of Pn-MOF. (A) Experimentally observed spectral evolution of the femtosecond TAS (excitation wavelength, 600 nm) of Pn-MOF in paraffin. 
(B) Temporal change of transient absorption at selected wavelengths and fitting curves from global analysis. (C) EAS obtained from global analysis based on the sequential 
model based on two components. EAS1 and EAS2 indicate the first and second components of EAS, respectively. (D) Temperature dependence of the transient absorption 


at 530 nm of Pn-MOF. 


within a hundred picoseconds (Fig. 3, A and B). Since the signal at 
530 nm is indicative of T|-T, absorption of pentacene chromophore 
(39, 45), we confirmed the ultrafast generation of the triplet pair 
state through SF. We can assign the second component as the triplet 
pair state with singlet multiplicity, "TT, because the spin multiplicity 
is conserved in the ultrafast timescale. We performed global analysis 
assuming a sequential model with two components to analyze time 
constants and spectral features (Fig. 3, B and C). The first and second 
components of the evolution-associated spectra (EAS) can reason- 
ably be assigned from their shapes to the spectra from Sı and 'TT, 
respectively. The S; state was converted to the 'TT state with a time 
constant of 16.8 + 0.1 ps, followed by negligible decay (>1 ns) in the 
current time window. The timescale of transition was slower than 
the typical timescale of SF in pentacene crystal (~100 fs) but still 
much quicker than the timescale of the typical intersystem crossing 
of pentacene in solution (>1 ns). The 'TT generation process likely 
outcompetes other photophysical processes such as intersystem 
crossing in a single chromophore. This is consistent with successful 
fine tuning of moderate intermolecular interaction required for 
efficient "TT generation. 

We then investigated the kinetics of the triplet pair state via 
nanosecond TAS in the Pn-MOF (Fig. 3D, red), which is crucial for 
QIS. Although the spectral shape was almost constant over time 
(fig. S5), multicomponent decay profile was observed. This complex 
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kinetics suggests that alternating spin multiplicity via interaction 
fluctuation affects the deactivation of the triplet pair state. To con- 
firm the origin of the multiple decays, we performed low-temperature 
(150 K) measurements (Fig. 3D, green), which clearly revealed 
slowed-down dynamics in the nanosecond timescale without altering 
the spectral shape (fig. S5), providing compelling evidence of the 
role of conformational fluctuation in the electron-spin coupled 
dynamics. 


Quintet multiexciton formation in Pn-MOF 

To quantitatively explore the spin dynamics in Pn-MOF, we per- 
formed CW-TREPR measurements of Pn-MOF in paraffin at room 
temperature. Spin polarization signals with typical A/A/E/A/E/E 
pattern of quintet (A and E represent microwave absorption and 
emission, respectively) and E/A pattern of triplet were observed on 
inner and outer side of the spectra, respectively (Fig. 4A). 

As mentioned above with the femtosecond TAS data, the EPR- 
silent 'TT state should be initially populated by SF. The spin conver- 
sion from ‘TT to “TT has been explained by the modulation of 
spin-spin exchange coupling J, where 6J = E('TT) — EČTT), with a 
conformational change in the covalently linked chromophore dimers. 
In this case, the chromophores remain strongly coupled to each 
other although the conformation of the dimer changes. In contrast, 
in dense organic crystals, triplet hopping forms weakly coupled 
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Fig. 4. Transient EPR spectra of Pn-MOF. (A) TREPR spectra for different delay times of Pn-MOF in paraffin at room temperature. Spectrum simulations by the rotation 
model analysis in Fig. 5 are shown by the red lines. (B) Conformations of the quintet TT; and TT>. The blue arrow represents dihedral angle (8) of the aromatic plane in Tg 
relative to T4 in the TaTg pair. (C) Decay profiles of the absorbance changes in nanosecond TAS (AOD) by the multiexcitons at temperatures of 300 K (red circles) and 150 K 
(blue circles), together with simulated profiles by the rotation model analysis (red and blue lines). Contributions of STT and SCTP (T +T) are shown by the solid and dashed 


lines, respectively, at 300 K (red line) and 150 K (blue line). 


triplet pairs, and mixing occurs between 'TT and °T'T, which are 
close to each other in energy level. The triplets then reencounter 
to produce a strongly coupled °TT. The MOF in this work is a unique 
system that combines the characteristics of both dimers and crystals 
because it not only has a structure in which many chromophores are 
accumulated but also has nanopores inside the crystal, which allows 
the chromophores to rotate and change their interactions. Therefore, 
it is appropriate to analyze the system using a model that assumes 
both J-modulation in the strongly coupled TT pair and spin mixing 
in the weakly coupled TT pair (Fig. 5). 

The °TT formation by the modulation of J through the chromo- 
phore rotations within the strongly coupled triplet pair in the 1-nm 
separated nearest pair shown in Fig. 2E (Fig. 5A) competes with the 
initial exciton migration to form the separated multiexcitons (Fig. 5, 
B and C). Depending on local exciton mobility by crystallinity or 
disorder, the weakly coupled triplet pair (Fig. 5B) and the non- 
correlated two triplets T + T (Fig. 5C) can separately be generated 
by the exciton migration. The quick triplet exciton dissociation 
at the crystalline domain would preferentially produce well-separated 
noncorrelated T + T and does not contribute to the electron spin 
polarization. In contrast, the weakly coupled triplet pair states are 
anticipated to be generated at the disorder region of the MOF to 
form the spin-correlated triplet pair (SCTP) (T + T) causing the 
nine adiabatic states (|1>, |2>, ---and |9>) by quantum superposi- 
tions of the nine of singlet, triplet, and quintet diabatic characters, as 
shown in Fig. 5B. These SCTP contributions are observed at the 
peak positions of the outermost magnetic fields (321 and 365 mT) 
exhibiting the E/A spin polarization pattern in Fig. 4A by the spin- 
spin exchange coupling of J = 1 mT in the 1.5-nm separated exciton 
pair (Fig. 2E). 
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From Fig. 5, we computed time development of the singlet, triplet, 
and quintet diabatic characters for the three different exciton-pair 
conformations [TT,, TT, and weakly coupled (T + T)] together 
with the different J values to reproduce the time-dependent EPR 
data using the density matrix formalism analysis (table $3). 

These diabatic spin characters in Fig. 5 are described by the den- 
sity matrix formalism of the equation of motions as follows: 


Prr, = — ier,» Prr, ]— (kiz + kpiss) Pr, (1) 
+ kpack Prat + ky Prr, + KrecPrr, 
Prr, = — ilHrr, Prr,] — ka Per, +Ki2Prr, (2) 


Drar = — Ary Prec] +kpissPrr, —kpackPrer —KrPrir (3) 


where p represents the density matrix for TT), TT, and T + T states 
and is numerically solved with the spin Hamiltonian f composed of 
the Zeeman interaction (rr), the zero-field splitting interactions 
(irns) in the individual two triplets, the spin-spin dipolar coupling 
(Hirss) between the triplets, and the exchange interaction (HrTex) 
as reported previously. Because H,-7,¢, and A;;,, are anisotropic and 
dependent on the TT conformations, the J values are also altered by 
the molecular motions, resulting in the anisotropic electron spin 
polarization in Fig. 4A caused by the conformation dynamics. This 
exchange process by kj) and kz; in Eqs. 1 and 2 is thus regarded as 
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Fig. 5. Model of the density matrix formalism analysis. Different multiexcitons, (A) strongly coupled triplet pair STT CTT; and 5TT)), (B) weakly coupled SCTP (T +T), 
and (C) noncorrelated two triplets T + T, are included with considering molecular rotation, exciton migration, and reencounter processes. 


the noise to fluctuate the J corresponding to the Debye exponential 
model (26). 

The EPR spectra were simulated as shown by the red lines in 
Fig. 4A. The electron spin polarization pattern (A/A/E/A/E/E) in the 
quintet is sensitive to the mutual orientations of the excitons in the 
TT and are generated by fluctuation of the exchange coupling. For 
the simulation of the EPR spectra, we started with the dihedral an- 
gle (P = 130°) between the aromatic planes in TT, as shown in 
Fig. 4B. We then carefully examined effects of the molecular confor- 
mation angles of Tg (eq. S14) in the T,-Tx pairs on the electron spin 
polarization of the quintet state. The dihedral angles (B) of the TT; 
or TT? states largely alter the spin polarization pattern (fig. S6). To 
examine whether the uniqueness of the motion is extracted from the 
spectrum, we also consider the angle distributions in the TT con- 
formations with fixing the TT, conformation (P = 130°) by averaging 
the five transverse magnetizations with weighting 1:4:6:4:1 ampli- 
tude for B = 130°, 140°, 150°, 160°, and 170° in the TT, as shown by 
the blue line in fig. S7A. In addition, averaging by the angle distribu- 
tion is performed in the TT; conformation, as shown in fig. S7B. Because 
the 130° to 150° spectrum (solid red lines in fig. S7) well agrees with 
the averaged spectra (dashed blue lines in fig. S7), it is concluded 
that the computation obtained by the 130° to 150° angle change may 
include the inhomogeneous distribution in the conformations with 
a large variation of +20°. Within this inhomogeneity, the average 
angle B = 150° is extracted for the Tg moiety in the TT, state as 
imposed in Fig. 4B, denoting that the conformation changes by mol- 
ecular rotations (red arrows in Fig. 4B) are coupled to the modulation 
of the exchange interaction in table S3. The stronger J (=—1.4 x 
10° MHz) in the TT, than that (J = —1.9 x 10 MHz) in the TT; is 
consistent with the dimer model average geometry of TT; and TT2 
causing contact edge-to-edge separations (0.29 and 0.30 nm, respec- 
tively) between carbon atoms in the Pn moieties in Fig. 4B based on 
the simulated structure (Fig. 2E). This supports the molecular 
rotation along with the Pn-ligating direction for the J-modulation. 
We would like to emphasize that we do not want to have a 
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quantitative discussion regarding the dihedral angle because of the 
heterogeneity in the MOF structure but rather to qualitatively un- 
derstand the importance of the suppressed motion with aver- 
age angles. 

Longer-lived multiexciton EPR signals are also obtained at 150 K 
(fig. S8). Experimental and simulated EPR spectra showed better 
agreement at 150 K than at 300 K. This may be due to the activation 
of several different low-frequency modes of conformational change 
at higher temperatures, making a perfect fit at 300 K difficult by 
considering only a single mode of conformational change. We per- 
formed the computation of the quintet EPR spectra with setting 
B = 125° and 0 = 40° in TT, with varying the TT2 conformation 
(fig. S9 and table S4). The center-field absorptive spectrum is 
obtained when we consider the motion between B = 125° in TT, 
and B = 110° in TT), as shown by the red line in fig. S9A. We also 
computed the time development of the transverse magnetization using 
the above conformation motion (fig. S10B) and added to the data in 
Fig. 4A to obtain the data (red lines) shown in fig. S10C. The better 
fit with the experimental results strongly supports that the highly 
heterogeneous conformation motions dominate the multiexciton 
dynamics. The time development of spin polarization coincides 
with the nanosecond TAS data, demonstrating that the “TT should 
be formed by the large and disorder conformational motions depicted 
in Fig. 4B. The reliability of the present numerical model analysis 
was confirmed by reproducing the EPR patterns of previous rigid 
tetracene dimers (fig. $11) (14). The EPR pattern of Pn-MOF could 
not be reproduced well without taking motility into account, con- 
firming that the pentacene moieties in the Pn-MOF have rotational 
dynamics. 

Asa result, the time evolution of the EPR spectra was successfully 
reproduced by the model of Fig. 5 in which the hopping disso- 
ciation (kpyss) and reencounter (kgack) to form the quintet TT, triplet 
3(T So), and singlet ITT via triplet-triplet annihilation of “Cl Ty 
are considered together with the J-modulation between TT, and 
TT, to generate the anisotropy in the quintet sublevel populations 
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(Fig. 4A, red lines). This analysis was also applied to understand 
the time development of the transient absorption signals in paraf- 
fin for f > 1 ns in the absence of the magnetic field. Notably, the 
profiles (solid lines in Fig. 4C) of the sum of the triplet excitons 
PTT: +°TT: +71 1S + TT, + TT: + SCTP (T + T) + noncorrelated 
T + T] are consistent with the deactivation kinetics of the spin- 
polarized exciton pairs with “TT and (T + T) in Fig. 4A. Quick 
nanosecond decays in Fig. 4C are originating from the singlet deacti- 
vation process (kgec = 5.6 X 10’ st at room temperature) of ITT that 
is in equilibrium with °TT and (T + T). The longer-lived transient 
absorption profile at 150 K than at 300 K is due to slower singlet 
deactivation (kee = 2.0 x 10’ s7}). It is possible to reproduce the 
nanosecond TAS decay profile even in a scenario that does not in- 
clude the J-modulation, where first (T + T) is generated by triplet 
diffusion and then “TT is generated by reencounter (fig. $12). How- 
ever, in this scenario, the "TT should still be hardly generated at 100 ns 
after photoexcitation, which contradicts the actual result that the EPR 
spectrum at 100 ns showed a strong °T'T-derived signal (Fig. 4A). This 
confirms that J-modulation plays a major role in the generation of °T'T. 


Room-temperature coherence control of quintet state 

To directly control the four spin qubits of the quintet multiexciton at 
room temperature, we applied a transient nutation microwave pulse 
sequence using an X-band pulsed EPR spectroscopy (Fig. 6, A and B). 
Although the quintet states are minor at the later delay time exceed- 
ing 1 ps in Fig. 4, A and C because of the quick deactivation in the 
strongly coupled TT pair to the ground state, the echo-detected EPR 
spectrum of the quintet is observed even at Taelay = 1.1 ps (Fig. 6C). 
This suggests that some of the strongly coupled TT states survive, 
while other dominant TT states undergo quick deactivation and 
dissociation by hot exciton characters with relatively large angular 
molecular rotation shown in Fig. 4B. These hot multiexcitons should 
also be subject to substantial transverse spin relaxation and not con- 
tribute to the echo signal. At the field strength of By = 348.1 mT 
(blue arrow in Fig. 6C), we obtained the microwave-induced spin 
coherence of the quintet as shown in Fig. 6D. This profile is succes- 
sively explained by the superpositions of the nutation frequencies 
(©ms) contributed by the quintet quantum gates of Ams = +1 in 
Fig. 6B (Fig. 5A) described by the following equation 


B 
Om, = : V/S(S +1) —m,(ms + 1) 


(4) 
where ms (= 0, +1, and +2) represent the quintet sublevel quantum 
numbers causing the transitions of Ams = +1 with S = 2 in Fig. 5A 
(detailed analysis method is described in table $5). A decoherence 
time constant of T2p = 150 ns is used. This is consistent with the 
transverse relaxation time of T, = 122 ns in fig. 16D. Bı is the micro- 
wave field strength in Eq. 1 and is set to be Bı = 0.45 mT. The nuta- 
tion frequency varies with the quantum gate: @ms = V6 guBB1/A for the 
gates (Qo,+1) between Qo and Q41 and ®ns = 2gugB1/A for the gates 
(Q+2,+1) between Q42 and Q41. Therefore, anisotropic quintet sub- 
level distribution by the S-Qms interconversion in Fig. 5A highly 
influences the frequencies and amplitudes of the quantum beats in- 
duced by the coherent microwave pulses at a field strength of Bo. The 
quantum beat profile is well reproduced by the dynamic model of 
the quintet populations considering subtle conformational motion 
between the dihedral angles of B = 130° and 122° in the TT, exci- 
ton pair for TT, and TT, respectively, which causes the fluctuation 
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in the through-space separations between 0.29 and 0.32 nm (fig. S10A). 
In the simulated profile (red lines in Fig. 6, D and E), a minor lower- 
frequency component (green line in Fig. 6D) is included by the gates 
between Q+2 and Q, , denoting that the two different nutation fre- 
quencies of 25 and 31 MHz are involved with opposite polarizations 
(Fig. 6B) in the experimental profile. 

The °T'T-derived signal observed in CW-TREPR was reproduced 
by a large dihedral motion between B = 130° and 150°, but such a 
large conformational change would cause a large fluctuation of the 
zero-field splitting interaction, leading to a severe quantum de- 
coherence (fig. $17). The pulse EPR nutation profile could not be 
reproduced by dihedral motion between B = 130° and 150° (green 
dashed line in Fig. 6E) but could be reproduced by the subtle confor- 
mational motion between the dihedral angles of p = 130° and 122°. 
This suggests that large and small conformational motions coexist 
as major and minor components, respectively, in the Pn-MOF. Since 
the Pn-MOF in this study was not very crystalline, we assume that 
majority of the MOF structure (source of CW-TREPR signal) was 
disordered, resulting in large conformational changes, while the con- 
formational change was suppressed in some minor domains with high 
crystallinity (source of pulse EPR echo signal). There are probably 
many disorders near the surface of the rounded nanoparticles (fig. S2), 
and the chromophores in this disordered region should cause the larger 
dihedral motion and the quicker dephasing in the quintet coherence 
and thus do not contribute to the spin echo signal. It is anticipated, 
however, that the ordered 22 cm™’ collective motions occur at the 
ordered domains most probably in the center core region of the MOF 
nanocrystals with preventions of the conformation motions. These 
minor TT pairs are possibly responsible for the weak echo signals. 

To corroborate this mode of motion, we additionally carried out 
terahertz spectroscopy of Pn-MOF (fig. S18). Although the frequen- 
cy of the large conformational change (14 cm7') was outside the 
measurement range, sharp Lorentzian lines were observed at around 
25 cm’ near the simulated frequency of the small conformational 
change (22 cm™'). Similar low-frequency motion modes have been 
observed in other MOFs in systems with restricted ligand mobility 
(46), supporting the idea that suppressive collective molecular mo- 
tion is occurring in the Pn-MOF. The sharp terahertz signals in 
fig. S18 indicate that the ordered collective modes (>22 cm”) originate 
from the motion of the entire framework via the metal-ligand coor- 
dination bonds. Although the detailed assignment of the collective 
mode is out of the scope of the current work, these phonon modes 
might couple with the local molecular rotations of the chromo- 
phores. For the singlet-quintet conversion in the strongly coupled TT 
with the negative exchange coupling, it is essential that the frequency 
of the J-modulation is closer to the frequency relevant to the 
quintet-singlet energy gap based on the Redfield theory, as reported 
previously (20, 26). In addition, the quintet-singlet energy gap is re- 
quired to be thermally accessible to get the quintet populations. In 
the present analysis, quintet-triplet energy gap is set to be ~1 THz 
corresponding to the —6J value (30 cm7!) that is close to the re- 
ported values (47, 48). In this case, the higher frequency vibrations 
(>50 cm!) in the terahertz spectrum should not contribute to minor 
quintet production because these frequencies are much higher than 
—6J. The broad and congested spectral component around 20 to 
30 cm™' should be responsible for the disordered motions causing 
the major echo-silent quintet EPR. Therefore, it is reasonable to hy- 
pothesize that the low-frequency Debye peak should contribute to the 
echo-silent quintet due to the disordered, nonstate-specific vibronic 


7 of 12 


PZOT ‘60 Yolk UO BIO‘sOUNTOS' mmm: sdyy woy popeoumoq 


SCIENCE ADVANCES | RESEARCH RESOURCE 


A Tdelay Techo = 200 ns Data acquisition BQ, 
Q,, ro 
Time Q, ~@& 
Laser x (ns) b (ns) Spin echo Q nd 
Q2 Q- 
c D Tdelay = 1.1 us 
= 
s 
2 
$ 
£ 
ee eo 
300 s20 s40 a60 380 400 j Nutation puise duain x (ns) i 
Magnetic field (mT) ` 
F 
1.5F T T aI Ta j Z 
40 —ß: 130 2122 ] 3 
3 we lee a. te E B: 130° 2 150° > | = 
5 °° gs | 8 
§ 00 £ ] 3z 
= Cc J m 
© -0.5 2 | I 
So : E- 
E -1.0 $ z 
= Ww 1 = 
-1.5 o. 
0O 100 200 300 400 0 5 10< 


Nutation pulse duration (ns) 


Tdelay (us) 


Fig. 6. Transient nutation of the quintet state at room temperature. (A) Pulse sequence of the transient nutation measurement. (B) Two main quantum gates detected 
by the pulsed microwave at By = 348.1 mT in the quintet sublevels. (C) Echo-detected EPR spectrum at Tgelay = 1.1 us. The durations of the first and second pulses are x = 8 ns 
and tz = 16 ns, respectively. Fitting result using eq. S15 and table S5 is shown as a red line. (D) Transient nutation profiles as a function of x in (A) at Tdelay = 1.1 ps using 
t2 = 16 ns. The simulated profile is shown by a red line considering the anisotropic sublevel populations. At the bottom, decomposition of the computed nutation data by 
the gates Qo,+1 and Q42,+1 are shown by purple and green lines, respectively. (E) Simulated profiles by the different rotation models. The red solid line corresponding to the 
red line in (D) is computed by the conformation motion between the dihedral angles of $ = 130° and 122° of the Tg molecule in the TaTg pair for TT; and TT, respectively. 
Green dashed line is obtained with considering the large dihedral angle change between $ = 130° and 150° in Fig. 4B. (F) Dependence of tuelay on the echo intensity (black 
line) employing x = 8 ns, t2 = 16 ns, and Techo = 200 ns with B4 = 0.45 mT together with delay time dependence of the CW-TREPR intensity at 351 mT in Fig. 4A (blue line). 
Fitting result using eq. S1 and table S5 is shown as a red line. 


coupling, while the ordered motion characterized by the Lorentzian 
peak at 25 cm™ should contribute to the minor long-lasting echo 
due to the ordered molecular motions with the suppressive confor- 
mation dynamics. While it is not the main topic and beyond the 
scope of the current work, we will test this hypothesis using MOFs 
with different crystallinity and structures. 

Specifically, the limited dihedral angle fluctuation in f is revealed to 
result in superposition of the transverse magnetization by the gates of 
Qo,ı (purple in Fig. 6, B and D) and Q-2,-1 (green in Fig. 6, B and D) 
at Bo = 348.1 mT from fig. S19A, while the magnetization by the Q+2,+1 
gates is minor for the larger angle-fluctuation case in fig. S19B. The hin- 
dered rotation thus leads to detections of the high- and low-frequency 
nutation components determined by œo = V6 SpByV/h and @_2=2gp1pBy/h. 
This is also consistent with the echo-detected EPR spectrum (red line 
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in Fig. 6C) obtained using this rotation model. These results imply 
that the quantum beating behavior can be altered by the degree of 
the structural fluctuation. 


DISCUSSION 

For quintet multiple excitons generated by photo-induced SF it has 
been unclear what kind of molecular motion is required to achieve 
both (i) effective quintet generation and (ii) noise suppressed qubits. 
Conformational motion is necessary for the spin conversion from 
ITT to °TT. However, if the change in orientation angle of the chro- 
mophore relative to the magnetic field is too large, then not only 
does exciton hopping cause dissociation and recombination between 
excitons but also “TT quantum decoherence occurs, making it 
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impossible to use as a qubit. On the other hand, low-mobility triplet 
pairs with small orientation angle changes prevent exciton hopping 
and recombination. In addition, the suppression of transverse relax- 
ation due to fluctuations in the zero-field splitting interaction makes 
it possible to slow down decoherence and function as qubits even at 
room temperature. 

To explicitly understand this trap characteristic for manipulation 
of the microwave-induced spin coherences, we plotted the echo in- 
tensity at 348.1 mT as a function of the delay time (Taelay) (Fig. 6F). 
The decay time (1.2 us) of the echo signal of the quintet is much 
longer than the decay time (0.2 us) of the quintet signals (blue line 
in Fig. 6F) detected by the CW-TREPR and TAS methods and is well 
fitted by the rotation model calculations (red lines in Fig. 6, C and F) 
using kgrec = 1.5 X 10° s~! at room temperature (table S5), demon- 
strating that the minor quintet species in some ordered domain are 
only responsible for the quantum coherence. This room-temperature 
°TT quantum coherence based on moderate suppression of chro- 
mophore motion in MOFs is in contrast to previous "TT generation 
based on exciton diffusion in dense molecular crystals, which is not 
immune to triplet separation and decoherence at room temperature. 
The usefulness of "TT as qubits has been demonstrated at cryogenic 
temperatures, but appropriate suppression of molecular motion in 
MOEFs leads to room-temperature drive of TT qubits. 

Previous studies treated that the quintet generations for the 
different sublevels are based on the Fermi’ golden rule (and by Red- 
field theory) by which the spin conversion is caused as the perturba- 
tion between the 'TT and °TT under the J-fluctuations (14, 49). This 
treatment is valid for the large positive J-coupling by which the “TT 
energy is notably lower and longer-lived with respect to the 'TT in 
the covalently linked chromophore dimer systems. In the present 
sample, however, we conclude that the j-coupling is negative and 
that its magnitude is smaller than the thermal energy. In this case, 
Fermi’s golden rule treatment is no longer valid because the 'TT is 
stable and thus interacting with the °TT during the measurements. 
This is well consistent with the quick quintet deactivation in nano- 
second timescales (Fig. 4C). This negative J is reasonably caused by 
the direct through-space orbital overlap between the pentacene 
chromophores in Pn-MOF structure (Fig. 4B). 

We incorporated the effect of the spin lattice relaxation within 
the quintet sublevels to reproduce the delay time dependence of the 
echo (red line in Fig. 6F) in eq. $17. The T; value of 3.8 us is longer 
than the deactivation time of the quintet polarization (blue line) 
in the TREPR measurement. This Tı (and T, = 150 ns) is another 
strong evidence for the ordered molecular motions dominating the 
very minor quintet component in the pulsed EPR detection. Fur- 
thermore, kgec = 1.5 ys! determined by the echo trace is much 
smaller than kre = 56 ps! obtained by the TREPR analysis. This 
slow recombination process well supports the suppressive confor- 
mation dynamics for restriction of the too much orbital overlap to 
lengthen the quintet lifetime in Fig. 6F. 

The high sensitivity of the pulsed EPR spectroscopy enables one 
to selectively observe the minor species when the phase memory 
time is long enough, while the nanosecond T> species originated 
from the major species are completely filtered out. It is also noted 
that the pulsed EPR spectroscopy filters out the short T} components 
because it is required to place a waiting delay time of at least 200 ns 
between the second microwave pulse and the first pulse (Fig. 6A) due 
to an unavoidable deadtime problem originating from the defense 
pulse introduced during the strong microwave irradiations, meaning 
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that T, of ~20 ns is out of the detection. The notable mismatches be- 
tween the CW-TREPR measurement (Fig. 4) and the pulsed EPR data 
(Fig. 6) at room temperature suggest the presence of more than one 
component in the MOF The one is a high-crystalline ordered compo- 
nent that gives a spin echo signal, and the other is a less-crystalline 
disordered component that does not give a spin echo signal and 
whose EPR signal mostly disappears within one microsecond. 

We revealed that the chromophore motion alters the phase and 
amplitude of the frequency component of its quantum beating. The 
control of multiple quantum gates is also extremely important for 
future developments in molecular quantum computing with multi- 
ple qubits, such as verification of CNOT gates. 

It is notable that not only the generation of quantum coherence by 
microwave irradiation has been demonstrated at room temperature 
by simulating chromophore motion in an MOF but also the details of 
molecular motion (change in molecular orientation) that enable the 
realization of room-temperature qubits have been found. In the pres- 
ent system, the low-mobility part that can generate °TT coherence is 
only a part of the MOE It will be possible to generate >IT qubits more 
efficiently in the future by searching for guest molecules that can in- 
duce more such suppressed motions and by developing suitable MOF 
structures. Precise control of the assembly structure and motion of 
the chromophores in MOFs is expected to lead to the development of 
innovative materials useful for the molecular-based QIS. 


MATERIALS AND METHODS 

All reagents were used as received unless otherwise noted. Dry 
tetrahydrofuran (THF) was prepared by treating with Molecular 
Sieves 4A 1/8, Wako. Liquid paraffin, n-butyl lithium in hexane, and 
zirconium tetrachloride (ZrCl,) were purchased from Sigma-Aldrich. 
p-toluene sulfonic acid, tin(II) chloride dihydrate (SnCl)-2H,0), 
sodium chlorite (NaClO2), sodium dihydrogen phosphate (NaH2PO,), 
and acetic acid were purchased from FUJIFILM Wako Pure Chemical. 
4-Bromobenzaldehyde, ethylene glycol, 6,13-pentacenedione, and 
2-methyl-2-butene were purchased from TCI. 


General characterization 

'H nuclear magnetic resonance (NMR) (400 MHz) spectra were re- 
corded on a JEOL JNM-ECZ400 spectrometer using tetramethylsilane 
(TMS) as the internal standard. Elemental analysis was carried out 
using a Yanaco CHN Corder MT-5 at the Elemental Analysis Cen- 
ter, Kyushu University. PXRD patterns were measured on a Bruker D2 
Phaser (Cu-Ka, 30 kV, 10 mA). For making the simulated structure 
and corresponding simulated PXRD pattern, Materials Studio software 
was used (Accelrys, Material Studio Release Notes, Release 4.2, Ac- 
celrys Software, San Diego, 2006). The geometry and unit cell were 
optimized by Materials Studio Forcite tool (50). Scanning electron 
microscopy image was collected by a Hitachi SU8000 with an acceler- 
ating voltage of 2.3 kV. Thermogravimetric analysis was conducted on 
a Rigaku Thermo Plus EVO2 under N2. N, adsorption isotherm mea- 
surements were carried out on a Bel BELSORP-max. UV-vis absorption 
spectra were obtained on a JASCO V-670 spectrophotometer. Fourier 
transform infrared spectra were recorded on a SHIMADZU IRTracer- 
100 with a Smiths DuraSamplIR II ATR device. 


Sample preparation 
PDBA was prepared (Fig. 7) according to our previously reported 
method (36). 
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Fig. 7. Synthetic route of PDBA. 


Synthesis of 2-(4-bromophenyl)-1,3-dioxolane (2) 

Sixty milliliters of a toluene solution of 4-bromobenzaldehyde (5.0 g, 
27 mmol), ethylene glycol (4.3 g, 54 mmol), and p-toluene sulfonic 
acid (0.17 g, 0.90 mmol) was refluxed with a Dean-Stark trap for 
8 hours. The solution was neutralized with 50 ml of KCO; aqueous 
solution and then extracted with ethyl acetate (3 x 50 ml). The 
organic layer was dehydrated by Na2SO,, filtrated, and dried under 
reduced pressure. The orange crude oil was purified using silica gel 
column chromatography (chloroform/hexane = 1:1), giving a 
compound 2 as yellow oil. Yield: 5.24 g (85%). 'H NMR (400 MHz, 
DMSO-d, TMS): 5 = 7.57 to 7.64 (d, 2H), 7.35 to 7.45 (d, 2H), 5.73 
(s, 1H), and 3.90 to 4.08 (m, 4H). 

Synthesis of 4,4'-(pentacene-6, 13-diyl)dibenzaldehyde (3) 
n-butyl lithium in hexane (2.5 M, 6.0 ml, 15 mmol) was dropped to 
a solution of 2 (3 g, 13.1 mmol) in dry THF (120 ml) at —78°C under 
N2 atmosphere. A 20-ml THF solution of 6,13-pentacenedione (1.24 g, 
4.0 mmol) was added after stirring for 1 hour at —78°C. The reaction 
was warmed to room temperature and stirred for 24 hours. A satu- 
rated aqueous NH,Cl (150 ml) was added to quench the reaction. 
The mixture was filtered and dried under vacuum. The dark-red oil- 
like crude was dissolved in 30 ml of THF and added SnCl2-2H20 
(4.2 g, 22 mmol) and 9 ml of concentrated HCI. To avoid the light 
exposure, the reaction mixture was covered by an aluminum foil 
and stirred at room temperature for 24 hours. The resulting solution 
was added 150 ml of H20 and extracted with chloroform (3 x 100 ml). 
The organic layer was dehydrated by Na2SO,, filtrated, and dried 
under vacuum. Silica gel column chromatography (chloroform/ 
hexane = 1:1) gave 3 as a dark purple powder. Yield: 1.61 g (82%). 
"H NMR (400 MHz, DMSO-dg, TMS): 5 = 10.32 (s, 2H), 8.32 (t, 4H), 
8.25(s, 4H), 7.90 to 7.95 (d, 4H), 7.84 to 7.88 (m, 4H), and 7.31 to 
7.35 (m, 4H). 

Synthesis of PDBA 

3 (1.61 g, 3.31 mmol) and 2-methyl-2-butene (35 ml, 331 mmol) were 
dissolved in THF (200 ml) at 0°C under Ar atmosphere. NaClO, 
(1.2 g, 13.2 mmol) and NaH2PO, (2.50 g, 20.9 mmol) in H20 (40 ml) 
were added to the suspension and stirred for 3 hours. The reaction 
was slowly warmed to room temperature and stirred for 24 hours 
under dark. After removing THF under reduced pressure, the residue 
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was diluted with water and washed with chloroform and methanol, 
giving compound PDBA as a dark purple powder. Yield: 1.05 g 
(65%). 'H NMR (400 MHz, DMSO-d,, TMS): 5 = 8.32 to 8.37 (d, 4H), 
8.25 (s, 4H), 7.85 to 7.88 (q, 4H), and 7.78 to 7.82 (d, 4H). Mass 
spectrometry (matrix-assisted laser desorption/ionization-time-of- 
flight mass spectrometry): mass/charge ratio [M+] calculated C36H2204, 
518.57; found 518.33. Elemental analysis for C3s5H2204-2H2O: calculated 
(%) H 4.73, C 77.97, and N 0.00; found (%) H 4.92, C 77.90, 
and N 0.13. 

Synthesis of Pn-MOF 

ZrCly (9 mg, 0.0385 mmol) and acetic acid (220 pl, 3.85 mmol) in 
2.26 ml of DMF were ultrasonically dissolved in a 6-ml screw vial. 
After stirring for 6 hours at room temperature, PDBA (20 mg, 
0.0385 mmol) was added, and the mixture was ultrasonicated again. 
The obtained purple solution was degassed two times by freeze- 
pump-thaw in a 10-ml ampule and flame-sealed under vacuum. The 
degassed sample was wrapped with aluminum foil to minimize light 
exposure and kept in an oven at 403 K for 48 hours. The mixture was 
cooled to room temperature over 72 hours, filtered, and soaked in 
methanol for 12 hours. After removing the solvent by filtration and 
drying for 5 hours at 323 K, dark green powder was obtained. 


TAS measurement 

The femtosecond TAS and nanosecond TAS experiments were per- 
formed using home-built pump-probe setups. A Ti:sapphire regen- 
erative amplifier (Spectra-Physics, Spitfire Ace) was used as the light 
source, (pulse duration, 120 fs; repetition rate, 1 kHz; pulse energy, 
4 mJ per pulse; central wavelength, 800 nm), seeded by a Ti:sapphire 
femtosecond mode-locked oscillator (Spectra-Physics, Tsunami) 
with a pulse duration of 120 fs, a repetition rate of 80 MHz, and a 
pulse energy of 10 nJ per pulse. For femtosecond TAS, one of the 
outputs from the amplifier was directed toward an optical parametric 
amplifier (Light Conversion, TOPAS), where the output wavelength 
was tuned to 600 nm for use as the pump pulse. The other output 
was focused on a sapphire crystal (3 mm in thickness), generating 
white light in the 450- to 750-nm range for use as the probe pulse. 
For nanosecond TAS, output pulses from an optical parametric 
oscillator (Light Conversion, NT242; central wavelength, 600 nm; 
pulse duration, 3 ns) using the third harmonic generation (THG) of 
a Nd:YAG laser as the excitation pulse were used as a pump pulse. 
The magic angle (~54.7°) was adopted between the pump and probe 
polarizations. The probe pulse that passed through the sample was 
dispersed by a polychromator (JASCO, CT-10; 300 grooves/500 nm) 
and detected by a multichannel detection system with a complemen- 
tary metal-oxide semiconductor sensor (UNISOKU, USP-PSMM-NP). 
The Pn-MOF sample was enclosed between quartz glass plates under 
an inert gas atmosphere. The TAS measurements were conducted in 
a semitransparent position on the Pn-MOF sample. 


CW-TREPR and pulsed EPR measurements 

TREPR measurements were carried out at room temperature using 
an X-band (9.682 GHz in the CW-TREPR and 9.589 GHz in the 
pulsed EPR) CW/fourier transform (FT) EPR spectrometer (ELEXSYS 
II E580) with a dielectric resonator using a quadrature detection. 
Samples were photoexcited by the second harmonics (532 nm) of a 
Nd:YAG laser (Continuum, Minilite II; full width at half maximum 
of ~5 ns). A laser depolarizer (SIGMA KOKI, DEQ 1 N) was placed 
between the laser exit and the microwave cavity. Repetition rate and pulse 
energy were 10 Hz, and 1 mJ per pulse, respectively. The CW-TREPR 
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spectra were obtained using a transient mode in the Bruker Elexsys 
spectrometer with a microwave power of 4.7 mW. The TREPR signals 
were directly amplified by a video amplifier using a control program 
(Xepr) for a digitizer of SpecJet-II. The TREPR data were smoothed 
for denoising using a binomial-weighted moving average to time course 
with a window of 300 ns for the spectra and the two-dimensional 
plot. The spin echo measurements in Fig. 6 were performed by using 
PatternJet-II controlled by the Xepr program. Pulse sequence of the 
Hahn echo detection (laser-Tgelay-™/2-Techo-M-Techo-echo) was used 
in the pulse-EPR measurements. The power of microwave pulse was 
set to maximize the intensity of the spin echo signal of oxidized pen- 
tacene species (fig. S16) observed in the absence of the laser irradia- 
tions following the CW-TREPR measurement. In this, the pulse 
durations were x = 16 ns and t, = 32 ns, respectively. The echo- 
detected field sweep measurements of the triplet pairs were performed 
by using Taelay = 1.1 HS, Techo = 200 ns, x = 8 ns, and f2 = 16 ns, and 
the magnetic field was varied. Transient nutation measurements 
were performed by varying x from 0 to 400 ns. 

Pn-MOF was put into 2-mm capillaries, with or without paraffin. 
The capillaries were degassed at room temperature with an oil pump 
and sealed with a flame. 

From Fig. 5, we computed time development of the singlet, triplet, 
and quintet characters for the different exciton-pair conformations 
together with the J values using coupled stochastic Liouville equa- 
tions, as reported previously (20, 51). In the details of the applied 
parameters for the multiexciton conformations, the kinetics are 
listed in tables S3 and S4 to compute the CW-TREPR data and the 
pulsed EPR data, respectively. The singlet, triplet, and quintet diago- 
nal terms in density matrix elements were used to obtain delay time 
dependences of the excitons (Fig. 4C) of piJ1, piv’, pt , plt2 


ss? Pss > Pingme? Pmgms? 
po pi4T, ptT, and Pme where m (= +1, 0, —1) represents the quan- 
tum number of the triplet sublevels. From Fig. 5, the following relations 


were used: “TT, = aP `TT, = BaPa “Epa Sp 


mmg mmg mm 
ITT =pyg TT =p rand SCTP =p" + DnP + Yi, PT 
(20, 51). 


Terahertz-TDS spectroscopy 

Terahertz spectra of MOF were measured by terahertz-time do- 
main spectroscopy (TDS) spectrometer (Advantest TAS7400) in a 
transmission setup. The frequency resolution of the spectra ob- 
tained by the system was 7.6 GHz. The number of accumulations 
was 1024 at room temperature. MOF in paraffin was measured with 
polyethylene cell, and solid MOF samples were measured without 
polyethylene cell. 
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